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Abstract. The effects of acoustic wave absorption, mode conversion and transmission by a sun- 
spot on the helioseismic inferences are widely discussed, but yet accounting for them has proved 
difficult for lack of a consistent framework within helioseismic modelling. Here, following a 
discussion of problems and issues that the near-surface magnetohydrodynamics hosts through a 
complex interplay of radiative transfer, measurement issues, and MHD wave processes, I present 
some possibilities entirely from observational analyses based on imaging spectropolarimetry. In 



particular, I present some results on wave evolution as a function of observation height and inclin- 
ation of magnetic field to the vertical, derived from a high-cadence imaging spectropolarimetric 
observation of a sunspot and its surroundings using the instrument IBIS (NSO/Sac Peak, USA). 
These observations were made in magnetically sensitive (Fe I 6173 A) and insensitive (Fe I 7090 
A) upper photospheric absorption lines. Wave travel time contributions from within the photo- 
spheric layers of a sunspot estimated here would then need to be removed from the inversion 
modelling procedure, that does not have the provision to account for them. 
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C\J ■ 1. Introduction 



Developments in sunspot seismology trace back to the original suggestion by[fh omas. Cram & N ye ( 1982) 



that interactions between sunspots and helioseismic p m odes could be used to probe the sub-surface structure 



of sunspots. The analyses that followed lThomas et alj (119821) focussed mainly on changes in the frequency 
wavenumber spectrum (v - k) and in the modal power distribution. These studies led to the discovery of 'ab- 
sorption' of p modes by sunspots dBraun. Duvall & Labontefl987l) : about 50% of the flux of acoustic wave 
energy impinging on a sunspot is not observed to ret urn to the quiet Sun. Development of several local heli- 
oseismic techniques, viz. the ring diagram analysis dml 1988b. helioseismic holography dLindsev & Braunl 
1990) ance helioseismology dDuvall. Jefferies. Harvey & Pomeranta ll993). has since brought 



in new ways of probing the subsurface structure and dynamics of sunspots. However, the question, viz. is a 
sunspot formed, in the sub-surface layers, of a monolithic flux tube or a cluster of flux tubes?, still remains 
to be answered. An answer to this question would also address the dynamics of heat and material flow in and 
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around sunspots, and hence would have far reaching implications for the magnetohydrodynamics of solar 
and stellar magnetism. 



Applications of time-distance helioseismology appeared as a promising avenue with its 3-dimensional 
tomographic images of flow and sound speed structures beneath sunspots dKosovichev. Duvall & Scherrei 
2000; IZhao. Kosovichev & Duvalll200ll) . These early results showed an increased sound speed region ex- 
tending from about 4 Mm down to about 18 Mm with a maximum change of about 1 - 2 %, while the near 
surface layers in the 1-3 Mm depth range show a decrease in sound speed of similar magnitude. The flow 
pattern dZhao et alJuOOll) consists of a shallow (1.5 - 3.0 Mm) converging flow that feeds a strong down- 
flow beneath the sunspot dDuvall. D'Silva. Jefferies. Harvey & Schoulll9 96) up to depths of about 5 Mm. 
Though these results have features indicative of the cluster model (Par kerll979h . new developments and im- 
provements in several different fr onts in local helioseismology have served to emphasise the inadequacy of 
such analyses dGizon et al.ll2010j). In contrast to results from time-distance helioseismology, studies based 
on phase sensitive holography (IBraun & LindsevluOOOl) have shown phase shifts of waves consistent with 
a faster propagation in the near surface layers, in direct correlation with the surface magnetic proxies (e.g. 
LOS magnetogram signals), and which decrease monotonically with depth becoming u ndetectable at layers 
deeper than about 5 Mm. Recent new ways of travel time measurements and inversions dMoradi et al Eoiot 
Svand a et al EoT 1 ) show that the moat outflows around sunspots extend much deeper (up to about 4-5 Mm). 
These new developments have brought to the fore the dominant direct interactions between acoustic waves 
and magnetic fields, which leave too large a signal in measurements to be treated with the conventional 
methods of seismic inversions that club such effects into thermal perurbations. 



The early contentions that p mode absorption of sunspots could be used to probe them, thus, have come 
around a full circle to the real i zation, through theore t ical attempts at explaining the above surf ace effects 
dCallv. Crouch & Braunll2003t l&ouch & Calivl |2005[ ISchunker, Braun. Callv & Lindsevll2006l) , that they 
are the very processes that need to be accounted for before we proceed further in the application of the later 
developed local helioseismic techniques. 



2. New developments: unreliability of old results 



The physical setting and nature of changes in the global structure and dynamics of the Sun is consistent 
with a conventional helioseismic analysis procedure, viz. a linear first order perturbation to the equilibrium 
structure of the Sun that the p mode frequencies effectively sense. However, such a treatment is much less 
adequate to probe the influence of sunspots in the near surface layers. The major sources of inadequacy in the 
local helioseismic analyses of sunspots, as have been gleaned from recent research, can be identified to arise 
from two basic causes: (i) inadequate understanding and modelling of the interactions between the acoustic 
waves and the sunspots, where magnetohydrodynamic effects dominate, and (ii) inadequate identification of 
the helioseismic observables due to complexities in the observation and measurement procedures themselves. 
However, much of what are known as 'surface magnetic effects' contain subtle inter-mixture of physical, 
measurement and analysis issues and hence have contributions from both the above causes. 
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2.1 Surface magnetic effects 



The neglect of direct magnetic effects due to the pressure and tension forces of the magnetic field on the 
wave speed, while inverting either the travel times or frequency shifts, is the foremost of issues arising from 
the cause (i) above. Interesting, but not yet fully understood, revelations on the seismic disguises of the dom- 
inantly near-surface interactions bet ween the magnetic field and acoustic waves came fo rth from analyses 
based o n phase-sensitive holography ( Braun & Li ndsev 20 00l): the 'showerglass effect' of[Lindsev & Braun 
(12005 ah and the 'inclined magnetic field effect' of lSchunker. Braun. Callv & Lindsevl (12005b . The former ef- 
fect refers to strong surface phase perturbations that the upcoming acoustic waves in active regions undergo 
resulting in impairment of their coherence, similar to the blurring of images seen through a commercial 
showerglass. These are measured as phases of the so called 'local control correlations' of ingressing (in- 
going) and egressing (outgoing) waves with wavefield observed at a particular point in active region (or 
sunspot) and are found to increase almost exponentially with magnetic field strength B. The 'inclined mag- 
netic field effect' pertains to the penumbral regions, where there are anomalous c hanges that depend o n the 
inclination angle of the magnetic field and the line o f sight angle of observations ( ISchunker et al 1 l2005l) . To 
correct for the showerglass effect ILindsev & Braunl d2005bl) devised a magnetic proxy, which is a complex 
amplitude that depends on B 2 and is nothing but the reciprocal of the appropriate local control correlation. 
'Corr ected measurements' fol low upon multiplying local surface signal with the above proxy. Such correc- 
tions dLindsev & Braunl2005bl) show that the sub-surface aco ustic anomalies disappear below a depth of as 5 
Mm, in contrast to the time-distance helioseismic inferences (IDuvall et al.l 1996; Ko sovichev et aL 2000). 



It is also likely that additional effects such as changes in the path length of the waves due to thermal 
expansion or contraction make significant contributions of either sign, and hence incorrect estimates of 
changes in sound speed from those in travel times. In particular, the path length changes associated with 
the Wilson depression and the propagating nature of (magneto-)acoustic waves (due to reduced cut-off fre- 
quency) would add contributions of opposite signs in the wave travel times. Clearly, neither thermal nor 
magnetic perturbations alone can explain the inferences. 



2.2 Observational issues: radiative transfer effects 



The altered thermal conditions in sunspots mean that the transfer of spectral line radiation is different from 
that in quiet Sun, and in the case of Zeeman sensitive lines the polarization and shape of the line interfere 
with the Doppler measurement procedure. With the added situation that the character of waves also are 
changed due to the magnetic field, radiative transfer effects manifest in Doppler velocity signals through 
subtle interaction of the above changes: the second basic cause [case (ii)] described earlier, viz. inadequate 
identification of the helioseismic observables due to complexities in the observation and measurement pro- 
cedures themselves. Here the helioseismic observable is the phase or travel time of a wave observed within 
a sunspot. Because the magnetic field lowers the acoustic cut-off frequency, and because it converts some 
of the incident acoustic waves into upward propagating ones confined to follow the field lines, the phases of 
waves measured within sunspots depend sensitively on the height within the line forming layers. Any line 
of sight angle dependent changes in (line) optical depth would then manifest as changes in the wave phases, 
i.e. different locations within a sunspot located at on off-disk center position would yield different phases for 
waves. This radiative transfer effect has been brought out clearly in an observational study of a sunspot in Ni 
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I line (6768 A) using the Advanced Stokes Polarimete r (ASP) at the Dunn Solar Telescope of the National 
Solar Observatory at Sac Peak, Sunspot, New Mexico ( Rajaguru et alj|2007 ). 



3. Helioseismic signatures of wave evolution in the observable layers: an example 

As an example of the near-surface effects discussed in the previous Section, I present here results of a 
study using a high cadence imaging spectropolarimetric observation of a sunspot and its surroundings in 
magnetically sensitive (Fe i 6173 A ) and insensitive ( Fe i 7090 A ) upper p hotospheric absorption lines. 



The results of this study have already been published (Rajag uru et al.l 12010) and we refer the readers to 
this original paper for a detailed account. We restrict ourselves here to a brief account of the observations 
and major physical implications of the results. The observations were made using the Interferometric BI- 
dimensional Spectrometer (IBIS) installed at the Dunn Solar Telescope of the National Solar Observatory, 
Sac Peak, New Mexico, USA. We observed a medium sized sunspot (NOAA AR10960, diameter as 18 Mm) 
located close to the disk center (S07W17) on 2007 June 8. Our observations involved scanning and imaging 
in all the Stokes profiles (/, Q, U, V) of magnetic Fe 1 6173.34 A and in Stokes / of non-magnetic Fe 1 7090.4 
A , with a cadence of 47.5 s. A 7 hr continuous observation was chosen for our analysis 

Line-of-sight (LOS) velocities of plasma motions within the line forming layers are derived from the 
Doppler shifts of line bisectors. We use 10 bisector levels with equal spacing in line intensity, ordered from 
the line core (level 0) to the wings (level 9), and derive 10 velocity data cubes, y,(x, y, f)(i = 0, 9), for each 
line. For the magnetic line, we us e the average of bisector velocities from the left (/ + V) and rig ht (/ - V) 



circular polarization (CP) profiles (San karasubramanian & Rimmele 20021 Idel Toro Iniestal l2003) and those 



from the / profile for the non-magnetic line. The 10 bisector levels span the height range within the line 
formation region in an unique one-to-one way. 



3.1 Instantaneous Wave Phases and Helioseismic Travel Times 

Instantaneous wave phases in the form of phase shifts 6<p, o(v)=Phase[V,(v)VQ(v)], where v is the cyclic 
frequency of a wave and V is the Fourier transform of u, due to wave progression between two heights 
corresponding to any one of the bisector levels i = 1,2, 9 and level (the top most layer) are calculated 



(Raja guru et al.ll2007l) . The 10 different data cubes from each line are run thro ugh a standard /?-mode time 



distance analysis procedure in center-annulus geometry (Raja guru et al.ll2004l) . Here, travel times for A = 
16.95 Mm are analysed, because, given the sizes of observed region (radius « 29 Mm) and the spot (radius 
« 9 Mm), this is the optimum A that facilitates distinguishing clearly the ingoing and outgoing waves in the 
sense of their interactions with the spot. Height dependent contributions to out- and ingoing phase travel 
times r + and r~ from within the line forming layers are determined using 6Tf = Tq - rf (i = 1, 9). 

We show in Figure 1 d<p&.o and St% , due to wave evolution within the region bounded by the wing (level 
8) and core (level 0) formation heights, against y. The v values marked in the panels of Figure 1 are the 
central frequencies of 1 mHz band filters used. Keeping in mind that (5</>8,o have contributions from a larger 
set of waves (as discussed above), results in Figure 1(a) for the magnetic line show a surprising amount of 
correlation between the two measurements, and moreover exhibit a strikingly similar y dependence. These 
results immediately reveal several interesting aspects of magnetic field - acoustic wave interactions: (1) 
first of all they confirm that helioseismic waves incident on the sunspot see themselves through to higher 
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a) Fe I 6173 A, CP Wings vs Core b) Fe I 7090 A, I Wings vs Core 
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Figure 1. Instantaneous phase shifts, <50 8 _ o (v), and changes in ingoing wave travel times, 6t^ q , due to wave propagation 
between the formation heights of wings (20 km) and core (270 km) of Fe i 6173 A (panel a), and of Fe i 7090 A {panel 
b) against y of B. 

layers of its atmosphere with a striking dependence on y: a coherent let through of incident waves happen, 
peaking around y « 30°, maintaining a smooth evolution of time-distance correlations; (2) remembering that 
CP profiles of the magnetic line have maximum sensitivities for velocities within vertical magnetic field, it 
is seen that a large fraction of waves propagating upward within such field are due to helioseismic waves 
originating at distant locations; and, (3) provide direct evidences that ingoing wave travel times would cause 
observing height dependent signals in flow inferences from travel time differences. 



a) (5T + ao , Fe I 6173 A vs Fe I 7090 A 




Figure 2. Changes in outgoing wave travel times, 6r^ Q , due to wave propagation between the formation heights of the 
wings (20 km) and the core (270 km) of Fe 1 6173 A (in red) and 7090 A (in green) lines as a function of y of B. {Panel 
b), a cartoon depicting wavefronts from acoustic sources beneath the umbra, wave paths and line formation heights (see 
text for details). 

Outgoing waves at a given measurement location, in general, would consist of those locally generated 
and those generated elsewhere undergoing reflection at the photosphere directly below it. These latter com- 
ponent would be seen in neither <50,-,o nor Srf Q , as they are evanescent at the observing height. For locally 
generated waves, circular wavefronts from a source, while their upward propagating parts see themselves 
up through the magnetic field, would cause outgoing wave correlations yielding distinct signatures in 6rf 
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(see Figure 2(b)). Results in Figure 2(a), for 8t% from both the magnetic and non-magnetic lines, do indeed 
provide such a diagnostic: outgoing waves starting at higher height (line core) within the sunspot atmosphere 
and reaching the quiet-Sun at the chosen A have shorter travel times than those starting at a lower height 
(line wings) and reaching the same quiet-Sun location; since this is simply not possible, the only explanation 
for this observation is the one contained in our previous sentence and illustrated in Figure 2(b), viz., outgo- 
ing wave time-distance correlations are predominantly due to waves directly from sources just beneath the 
sunspot photosphere when oscillations observed within it are used. 



4. Discussions and conclusion 



Almost all time-distance helioseismic analyses proceed under the working assumption that wave signals at 
observation heights are evanescent and hence oppositely directed wave paths involving photospheric reflec- 
tions at two separated points are of identical path length. This assumption is basic to the inferences on flows 
and wave speed from travel time differences and mean, resp ectively . In an e arly theoretical study, accompan- 
ied by attempts to model the helioseismic observations of iBraunl ( 1997 ). Bogdan. Braun. Lites & Thomas! 
(1998) showed the influences of both the p-mode forcing of, and spontaneous emissions by, sunspots on 
acoustic wave travel times. Our analyses here have yielded transparent observational proofs for both effects, 
for the first time, with important new perspectives: (1) the process of transformation of incident acous- 
tic waves into propagating (magneto)-acoustic waves up through the magnetic field happen in a coherent 
manner allowing a smooth evolu tion of time-distance correlations and, in agreement with s everal recent 
theoretical and numeric al studies (Cally 200l ICrouch & Callvll2005l: ISchunker & Callvll2006l) . this process 
depends on the inclination angle (y) of magnetic field to the vertical, and (2) outgoing waves from acoustic 
sources located just beneath the sunspot photosphere add important additional contributions for both mean 
travel times and differences. Our results have also shown observational prospects for consistently account- 
ing for the above effects in sunspot seismology, viz. the indispensability of imaging spectroscopy to extract 
wave fields so as to be able to correctly account for the wave evolution within the directly observable layers 
of sunspot atmosphere. 
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